Cooling process for inelastic Boltzmann equations 
for hard spheres, Part I: The Cauchy problem 

S. Mischler\ C. Mouhot^ and M. Rodriguez Ricard^ 

February 2, 2008 

Abstract 

We develop the Cauchy theory of the spatiahy homogeneous inelastic 
Boltzmann equation for hard spheres, for a general form of collision rate which 
includes in particular variable restitution coefficients depending on the kinetic 
energy and the relative velocity as well as the sticky particles model. We prove 
(local in time) non-concentration estimates in Orlicz spaces, from which we 
deduce weak stability and existence theorem. Strong stability together with 
uniqueness and instantaneous appearance of exponential moments are proved 
under additional smoothness assumption on the initial datum, for a restricted 
class of collision rates. Concerning the long-time behaviour, we give conditions 
for the cooling process to occur or not in finite time. 
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1 Introduction and main results 

In this paper we address the Cauchy problem for the spatially homogeneous Boltz- 
mann equation modelling the dynamic of a homogeneous system of inelastic hard 
spheres which interact only through binary collisions. More precisely, describing the 
gas by the probability density f{t, f ) > of particles with velocity v G (A^ > 2) 
at time t > 0, we study the existence, uniqueness and the qualitative behaviour of 
solutions to the Boltzmann equation for inelastic collision 

(1.1) ^ = QifJ) in (0,+oo)xM^, 

(1.2) /(O,-) = /in in M^. 

The use of Boltzmann inelastic hard spheres-like models to describe dilute, rapid 
flows of granular media started with the seminal physics paper ^j, and a huge 
physics litterature has developed in the last twenty years. The study of granular 
systems in such regime is motivated by their unexpected physical behavior (with the 
phenomena of collapse -or cooling effecf- at the kinetic level and clustering at the 
hydrodynamical level), their use to derive hydrodynamical equations for granular 
fluids, and their applications. 

From the mathematical viewpoint, works on the Cauchy problem for Boltzmann 
models have been flrst restricted to the so-called inelastic Maxwell molecules model, 
where, in a similar way to the Maxwell model in the elastic framework, the colli- 
sion rate is independent on the relative velocity. Existence, uniqueness of solutions 
and description of the asymptotic cooling has been obtained in [7] for the inelas- 
tic Maxwell molecules model with constant normal restitution coefficients as well as 
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with some cases of normal restitution coefficients depending on the kinetic energy of 
the solution. More precise properties of the solutions, such as their convergence to 
self-similarity, have also been investigated and we refer to our companion paper |2Z| 
for more details and more references on this issue. At least in the spatially homo- 
geneous setting, the inelastic Maxwell molecules model seems well understood now. 
The Maxwell molecules model is important because of its analytic simplifications 
(with regards to the hard sphere model) allowing to use powerful Fourier trans- 
form tools as introduced by Bobylev (see for instance [S]) for the elastic Maxwell 
molecules Boltzmann equation. Another simplification which has lead to interesting 
results is the restriction to one-dimensional models (in space and velocity), where, 
on the contrary to the elastic case, the collision operator has a non-trivial outcome. 
These models have been considered in for some cases of normal restitution 

coefficients possibly depending on the relative velocity. 

It is possible to modify the collision operator of the inelastic Maxwell molecules 
model by a multiplication by a function of the kinetic energy in order to restore 
its dimensional homogeneity (see [7 ) and thus the rate of cooling. Physically the 
derivation of this model amounts to replace the collision rate by a mean value 
independent on the relative velocity, starting from the inelastic hard spheres model, 
and the resulting approximation is named pseudo-Maxwell molecules in [7] . However, 
fine properties of the distribution (such as the behavior of the overpopulated tails 
or the self-similar solutions) are broken or modified by that approximation with 
respect to the inelastic hard spheres model. The recent papers [18t |S| have studied 
the case of inelastic hard spheres with constant normal restitution coefficients in any 
dimension and in various regimes: in particular in a thermal bath, i.e., when a heat 
source term is added to the equation, and in the self-similar variables of the free 
regime. Existence and smoothness of solutions to the Cauchy problem and to the 
associated stationary problem are obtained in ^H] for the thermal bath regime, while 
precise estimates on the tails of the stationary solutions (assuming their existence) 
for various regimes (including the two ones above-mentioned) are exhibited in [S]. 

In the present work, we shall construct solutions to the freely cooling Boltz- 
mann equation for inelastic hard spheres in any dimension N > 2 and for a general 
framework of distributions of inelasticity (defined by a measure on the set of all 
possible post-collisional velocities), which covers in particular variable normal resti- 
tution coefficients possibly depending on the relative velocity and the kinetic energy 
of the solution. It includes the cases of visco-elastic hard spheres model (see jH]) as 
well as the case of sticky particles model. Our framework enables to consider inter- 
esting physical features, such as elasticity increasing when the relative velocity or 
the temperature decrease ( "normal" granular media) or the opposite phenomenon 
( "anomalous" granular media). We refer to jTHHElElini and the references therein 
for a physical motivation. Let us emphasize that these solutions are new even in 
the case of a constant normal restitution coefficient as considered in |18| IS]. We also 
discuss the uniqueness of solutions, the instantaneous appearance of exponential mo- 
ments and various conditions on the collisions rate for the collapse to occur or not 
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in finite time. A second part of tliis work [27^ will be concerned with the existence 
of self-similar solutions and the tail behavior of the distribution. In a third part 
[2H1, we shall prove the uniqueness and the asymptotic stability of these self-similar 
solutions for a small inelasticity. 

From the viewpoint of mathematical tools, our main new contributions can be 
summarized as follows: 

(i) A generalization of the propagation of the L^-norm of the solution for the 
elastic Boltzmann equation based on Young's inequality as introduced in ^2] (see 
also IHOl QE| where similar ideas are used for a different model), into a result of 
propagation of Orlicz norms for inelastic (and elastic) Boltzmann equations. This 
a priori estimate is used in order to prove the existence of solutions to the inelastic 
Boltzmann equation with energy dependent inelasticity. Let us emphasize that it 
also gives an alternative proof of existence of solution for the elastic Boltzmann 
equation with initial datum having only finite mass and kinetic energy (but possibly 
infinite entropy). 

(ii) A generalization of the DiBlasio uniqueness Theorem for the elastic hard 
spheres Boltzmann equation (see jTUl 1201 IHOl IM] ) and for the inelastic hard spheres 
Boltzmann equation with constant normal restitution coefficients (see jlSlEI), to 
the inelastic hard spheres Boltzmann equation with energy dependent normal resti- 
tution coefficients (see also [121 where similar tools are developed). 

For points (i) and (ii), one of the main ideas of the proof is an appropriate use of 
the change of variables v' (for fixed {v,a)) and v ^ v' (for fixed (f*,(j)) in 

the spirit of the proof of the so-called "cancelation lemma" introduced in [HTj (see 
also |1J). 

(iii) An improvement of the result of propagation of exponential moments for the 
elastic Boltzmann equation jB] and for the inelastic Boltzmann equation |S], into a 
result of instantaneous appearance of exponential moments. This is obtained by 
combining estimates from [S] together with a simple o.d.e. argument introduced in 
the context of the Boltzmann equation in |4Uj . 

Before we explain our results and methods in details, let us introduce the prob- 
lem. 



1.1 A general framework for the collision operator 

We denote by B the rate of occurance of collisions of two particles with pre-collisional 
velocities which gives rise to post-collisional velocities The collision 

may be schematically written 



;i.3) W + K}^M + K} with 



v' + v' = V + V:^, 
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More precisely, for any fixed pre-collisional velocities v, G M^, we introduce a 
parametrization by z G -D := {w G M^; \w\ < 1} of all possible resulting post- 
collisional velocities in p.3p in the following way: 

{v' = {v + v^)/2 + z\v^-v\/2 
< = {v + v^)/2- z\v^-v\/2. 
The collision rate B takes the form 

(1.5) B=\u\h, b = a{£)P{£,u;dz) 

where u = v — is the relative velocity, a is an intensity coefficient, P is the 
normalized cross-section (it is a probability measure on D for any fixed S,u), and 
£ is the kinetic energy of the distribution /, defined by 

£ :=£{/)= [ f\v\'dv. 

JRN 

The expression ()1.5|) refiects the fact that we are dealing with hard spheres which un- 
dergo contact interactions. The term \u\ C({£) corresponds to the rate of collisions of 
two particles with pre-collisional velocities v, G M^, while the term j3 corresponds 
to the conditional distributional probability to obtain the two post-collisional veloc- 
ities {v',v'^}. The non-negative real 1^1 is the restitution coejficient which measures 
the loss of energy in the collision, since 

(1.6) \v'\^ + \v'f - \v\^ - \v,\^ = -^(1 - \z\^) \v, - < 0. 

In the above formula, |z| = 1 corresponds to an elastic collision while z = corre- 
sponds to a completely inelastic collision (or sticky collision). 

The bilinear collision operator Q{f, /) models the interaction of particles by 
means of inelastic binary collisions (preserving mass and total momentum but dissi- 
pating kinetic energy). More precisely, we define the collision operator by its action 
on test functions (which is related to the evolution of the observables of the probabil- 
ity density). Taking ip = ip{v) to be some well-suited regular function, we introduce 
the following weak formulation of the collision operator (valid under the symmetry 
assumption ()1.11|) below on /3) 

(1.7) (Q(/,/),^) -.= 1 ! [ hf ! {^: + ^' -^-v,)B{£,u;dz)dvdv,. 

Here and below we use the shorthand notations ip '■= i^i^), ip* '■= "^("y*), i^' '■= i^iy') 
and Tp'^ := ip{v'^) for any function ip on M^. 

A first simple consequence of the definition of the operator ()1.7p and of the 
parametrization p.4|) is that mass and momentum are conserved 

dtj^N-'yvj 
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a fact that we easily derive (at least formally), multiplying the equation (jl.lj) by 
if = 1 OT if = V and integrating in the velocity variable (using (jl.7j) ). In the same 
way, multiplying equation (jl.ll) by ip = |f p, integrating and using ()1.6|) and ()1.7p . 
we obtain that the kinetic energy is dissipated 

(1-8) j/{t) = -D{f) < 0, 

where we define the energy dissipation functional D and the energy dissipation rate 
A, which measures the (averaged) inelasticity of collisions, by 

u) dvdv^, 



D{f) ■■= [ [ fMufA{S, 

A{£,u) := - I {l-\zf)b{S,u;dz)>0. 
4 Jd 



Finally, we introduce the cooling time, associated to the process of cooling (pos- 
sibly in finite time) of granular gases: 

(1.9) Tc := inf |t > 0, S{t) = Vt > t| = sup > 0, S{t) > Vt < sj. 

This cooling effect (or collapse) is one of the main motivations for the physical and 
mathematical study of granular media. 

The Boltzmann equation (ll.lj) is complemented with an initial condition ()1.2|) 
where the initial datum is supposed to satisfy the moment conditions 

(1.10) < /in G L;(M^), / findv = l, [ fi^vdv = 

for some q > 2. Notice that we can assume without loss of generality the two last 
moment conditions in ()1.10|) . since we may always reduce to that case by a scalling 
and translation argument. Here we denote, for any integer g e N, the Banach space 

Ll = \f:R^ — > M measurable; ||/||ii := /" \f{v)\{l + \v\'')dv<oc\. 
We also define the weighted Sobolev spaces W^'^ (g G M and A; G N) by the norm 

\s\<k 

We introduce the space of normalized probability measures on M^, denoted by 
M^(M^), and the space BVg{R^) (g G M) of (weighted) Bounded Variation func- 
tions, defined as the set of the weak limits in of sequences of smooth functions 
which are bounded in W^'^iM.'^). Throughout the paper we denote by "C" various 
constants which do not depend on the collision rate B. 
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1.2 Mathematical assumptions on the coUision rate 

Let us state the basic assumptions on the coUision rate B: 

• The probabihty measure (5 satisfies the symmetry property 

(1.11) P{£,u;dz) = (3{£,-u;-dz). 

• For any 9? G Cc(M^) the functions 

(1.12) {v,v^,£)i-^ / (p{v') f3{£,u;dz) and £ i-^ a{£) 

Jd 

are continuous on M.^ x M.^ x (0, 00) and (0, 00) respectively. 

• The probabihty measure (3 satisfies the following angular spreading property: 
for any £^ > 0, there is a function js{^) > such that 

(1.13) V£ > 0, e / f3{£, u; dz) < jeie) 

J{\u-z\e[-l,l]\[-l+s;l-e]} 

and jsis) ^ as e — uniformly according to £ when it is restricted to a 
compact set of (0, +00). 

We will sometimes restrict our analysis to a kind of generalized (energy depen- 
dent) visco- elastic model assuming that the cross-section b reduces to an absolutely 
continuous measure according to the Hausdorff measure on the sphere 

(1.14) = 

More precisely, we assume that 

(1.15) b{£, u; dz) = 5{2=(i-e)n/2+(i+eW2} b{£ , \u\,u- cr) da 

where da is the uniform measure on the unit sphere, 6 is a non-negative measurable 
function and e : (0, 00) x x [—1, 1] [0, 1], e = e{£, \u\,u ■ cr) is a continuous 
function. For a vector x G ]R^\{0}, we define x = x/\x\ and §^^^ stands for 
the unit sphere of M.^ . Roughly speaking, the generalized energy dependent visco- 
elastic model corresponds then to the case where for any direction z G E>^~^, the 
post-coUisional velocities {v',v'^) such that {v' — vl)/\v' — v'J — z are uniquely (or 
deterministically) defined by the pre-coUisional velocities (v,v*). 

For the uniqueness of the energy coupled models, we shall need the following 
additional assumption: 
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HI. The cross-section b satisfies (jl.l5j) with b bounded, e = e{S) and the following 
locally Lipschitz conditions holds: for any compact subset K C (0, oo) there 
exists a constant Lk G (0, oo) such that for any S,S' & K 

(1.16) sup \\b{S', u, .) - b{S, u, .) ||Li(§iv-i) <Lk\S'-S\ 
and 

(1.17) \e{S')-e{S)\<LK\S'-S\. 

In the study of the cooling process, we always assume: 

H2. The energy dissipation rate A{S,u) in ()1.9|) is continuous on (0, +oo) x 
and satisfies 

(1.18) A{£,u)>0 yueR^,£>0. 

We will also need one of the two following additional assumptions: 
H3. For any £o,Soo ^ (0, oo) (with Sq > S^o) there exists ip such that 

(1.19) A{S,u)>^{\u\) G (^oo,^o), Vm G R^, 

with ip G C(R_|_, ]R+) and such that for any R> there exists ipj^ > with 

(1.20) ipi\u\) >iPr\u\~^ \/u eR^, \u\> R/2. 

This assumption is quite natural. In particular, it holds for a "normal" gran- 
ular media. 

H4. The cross-section b satisfies ()1.15|1 with e = e{S, \u\) and there exists bo,bi G 
(0, oo) such that bo < b < bi a.e. and x ^— b{S, \u\,x) is nondecreasing and 
convex on (—1, 1) for any fixed S G (0, oo) and u G M^. 

Notice that under assumption (jl.lSp with b = b{u ■ a) and e = e{S, u) the energy 
dissipation rate just writes 

(1.21) A(£:,M) = C^(l-e2), 

where Cat is a constant depending on the dimension. 

Let us emphasize that the classical Boltzmann collision operator for inelastic 
hard spheres with a constant normal restitution coefficient e G [0,1], as studied 
in j7] and jTHj, is included as a particular case of our model, and satisfies all the 
assumptions above. But the formalism described from p.3p to p.l3p is much more 
general than this case. In particular, we may also consider: 

1. Uniformly inelastic collision processes such that 

(1.22) 3^oe[0,l) s.t. supp5(^,u,.) C ^(0,20) VnGR^, V^>0, 
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which includes the sticky particles model when zq = 0. 

2. The physically important case ()1.14pi.l5|l of collisions defined by a normal 
restitution coefficient e and the cross-section b which possibly depend on 8, \u\ and 
u ■ (7. In particular it covers the kind of models studied in [7j (where e depends 
on S, and b is independent on S and \u\). It includes also the important case of 
the visco-elastic hard spheres model where b = b{u ■ a) and the normal restitution 
coefficient depends (smoothly) on the normal component of the relative velocity, 
that is 1^1 Im — a\/2 in our notation (see P). 

3. This formalism also covers multidimensional versions of the kind of models 
proposed in , which corresponds to the case where b is the product of a measure 
depending on 1^1 and a measure of d ■ z absolutely continuous according to the 
Hausdorff measure. One easily checks that our assumptions p. 511. 1111. 1211.1^ on 
the collision rate are quite natural for this kind of models as well. Note that our 
measure framework for B can also models situations where, in the opposite to the 
generalized visco-elastic case, there is some stochasticity or uncertainty on the degree 
of inelasticity of the collisions, for instance due to some experimental noise, or due 
to the fact that particles in the g mixture of different inelasticity behaviors, 
which are therefore handled statistically. 



The fact that 6 is a ffnite measure on D allows to deffne the splitting Q 
where and Q~ are defined in weak form by 



(1.23) 

and 
(1.24) 



f if' \u\b{£,u;dz) dvdv^, 
Jd 



{Q {9,f),v) ■= I I 9*f I v\u\b{8,u]dz)dvdv^, 



where v' is defined by (II. 4|) . A straightforward computation shows that it is possible 
to give a very simple strong form of Q~ as follows 

(1-25) Q-igJ) = Lig)f, 

where L is the convolution operator 



(1.26) 



L{g){v) := a{£) / g{v^) \v - dv^ 



Under assumption ()1.15|) . the expression of reduces to 
(1.27) {Q'^ig, f),y^) ■= fl'*/l'"l / if'b{S,\u\,u-cr)dadvdv^ 

where v' is defined by the formula (deduced from ()1.4|1 and ()1.14|l ) 



1.28) 



V = V — 



1 + e 



u 



\u\a 



1 + e 



u 



\u\a 
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1.3 Statement of the main results 

Let us now define the notion of solutions we deal with in this paper. 

Definition 1.1 Consider an initial datum /in satisfying M.10\} with q = 2. A non- 
negative function f on [0,T] xR^ is said to be a solution to the Boltzmann equation 
dll-dll if 

(1.29) /gC([0,T];L1(M^)), 

and if ^H^-^^ holds in the sense of distributions, that is, 

(1-30) rU f^dv{Q{fJ),<P)]dt=[ f,nm-)dv 

for any(f)eCl{[0,T) x M^). 

It is worth mentioning that (jl.29|) ensures that the collision term Q{f,f) is 
well defined as a function of L^(]R^). Indeed, on the one hand, we deduce from 
/ G C([0,r];Li(R^)) that £{t) G Ki on [0, T] and thus a{£{t)) G K2 on [0,T] for 
some compact sets C (0,oo). On the other hand, from the dual form ()1.23|) it 
is immediate that is bounded from L\ x L\ into L^, with bound a{S) (see also 
fHlEZ! for some strong forms of the Q~^{f, f) term). It turns out that a solution /, 
defined as above, is also a solution of p.l|l - ()1.2j) in the mild sense: 

f{t,-) = fin+[ Q{f{s,-),f{s,-))ds a.e. in M^. 
Jo 

Another straightforward consequence is that if / G L°°{[0,T), L^) then / satisfies 
the following chain rule 

(1.31) ^[ E{f)<pdv = {Q{f,f),E'{f)<p) in V'{[0,T)), 

at j^N 

for any S G C\R) n iy^'~(M), G L^q(M^), in the sense of distribution on [0, T). 

Let us state the main results of this paper. First, we give a Cauchy Theorem 
valid when the collision rate B is independent on the kinetic energy. 

Theorem 1.2 Assume that B satisfies the assumptions M.^) - n~Tl\} - m. l'^) - M. 1 5j) 

with h = h{u] dz): the cross-section does not depend on the kinetic energy. Take an 
initial datum /in satisfying M.IO^) with q = 3. Then 

(i) For all T > 0, there exists a unique solution f G C{[0,T];LD fl L°°{0,T; L^) 
to the Boltzmann equation l[l.l]} - n~^} . This solution conserves mass and mo- 
mentum, 

(1.32) / f{t,v)dv=l, / f{t,v)vdv = VtG[0,T], 

and has a positive and decreasing kinetic energy 

(1.33) < S{t2) < £{ti) < fin = f (0) Vti G [0, T], ti< t2. 

In particular, the life time of the solution (as introduced in / li.^^) ) is T^ = +00. 



10 



(ii) Moreover, assuming H2-H3 or H2-H4 (with e and b independent on the 
kinetic energy), there holds 

(1.34) S{t)^0 and /(t, .) ^ 6^=0 in M^{R^)-weak* when t ^ T^. 

In other words, the cooling process does not occur in finite time, but asymp- 
totically in large time. 

Remarks 1.3 Let us discuss the assumptions and conclusions of this theorem. 

1. Under assumption H4 and when the collision rate is independent on the ki- 
netic energy, one can prove in fact that there exists a unique solution f G C([0, 00); L^) 
satisfying and for any initial condition /in satisfying il.lU^) with q = 2. 
The proof is quite more technical and we refer to ^31 J where the result is presented for 
the true elastic collision Boltzmann equation; nevertheless the proof may be readily 
adapted to the inelastic collisional framework. 

2. The existence and uniqueness part of Theorem \1.2\ (point (i)) extends to a 
collision rate B = B{u; dz) > which satisfies the sole assumptions 



for some constant Co G This corresponds to the so-called cut-off hard potentials 
(or variable hard spheres) model in the context of inelastic gases. 

3. For a uniformly dissipative collision model, i.e., such that 



a fact which holds under assumption or under assumption H4 with a normal 

restitution coefficient e satisfying e{\u\) < Cq G [0, 1) for any u G M^, we may prove 
the additionnal a priori bound 



As a consequence, one can easily adapt the proof of existence and uniqueness in 
Theorem \1.2i and then one can easily establish that the existence part of Theorem \1.2\ 
holds for any initial datum /in satisfying M.1U\) with q = 2. 

4. The existence and uniqueness part of Theorem \1.2\ (point (i)) immediately 
extends to a time dependent collision rate B = \u\ 7(t) b{t, u; dz) where b{t, u; ■) is a 
probability measure for any u G MJ^ , t G [0, T] such that b{t, u; dz) = b{t, —u; —dz), 
and '-/{t) is a non-negative function in L°°{0,T). 




A{u) > Ao G (0,00) 
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5. Finally let us emphasize that Theorem li.jjl applies to the important ( non- 
coupled) model of visco- elastic hard spheres. Indeed the collision rate of this model 
satisfies assumptions M .5M . 1 111 . . 1'J^) as well as H2 and H3, with b and e in- 
dependent of S. We refer to the work in preparation 129^ which shall be devoted to 
the detailed study of this particular model. 

Now, let us turn to the case where the colhsion rate depends on the kinetic 
energy of the solution. 

Theorem 1.4 Assume now that B satisfies the assumptions M.^) - 1[ni\] - M.l'^) - 
M.13\) and that the cross-section b = b{£, u; dz) indeed depends on the kinetic energy 
£. Take an initial datum /in satisfying M.1U\] with q = 3. 

(i) There exists at least one maximal solution f G C([0,T];L2) fl L'^iQ^T] L\) , 
VT G (0,Tc), for some Tc G (0, +cx3], which satisfies the conservation laws 
and the decay of the kinetic energy 

(a) If the collision rate satisfies the additional assumption HI, and the initial 
datum satisfies the additional assumption /in G BV4 fl L^, then this solution 
is unique among the class of functions C{[0,T], L^) fl L^{0,T; Ll) , for any 
TG (0,T,). 



(Hi) The asymptotic convergence ^1.^4 ^ holds under the additional assumptions 
H2-H3 or H2-H4. 

(iv) If one of following assumptions a. or b. is satsfied, then Tc = +00; 

a. a is bounded near 8 = and js converges to as e ^ uniformly near 
E = 0; 

b. B satifies H4, A is bounded by an increasing function Aq which only 
depends on the energy, and /in e''" 1^1" G with r] G (1,2], > 0. 

(v) IfA{£,u) > AqS^ with Ao>0 and 6 < -1/2, then T^ < +00. 

Remark 1.5 Under the assumptions of point (ii) on the initial datum, by using a 
bootstrap a posteriori argument as introduced in ]31^. one can prove that there exists 
a unique solution f G C([0, cxd); L^) satisfying M.H^) and for any initial 

condition /in satisfying M.l(^) with g > 4 and /in G BV^. 



1.4 Plan of the paper 

We gather in Section |21 some new integrability estimates on the collision operator 
which can be of independent interest. We prove convolution-like estimates in Or- 
licz spaces for the gain term. We give then estimates on the global operator in 
Orlicz space, which show essentially that even if the bilinear collision operator is 
not bounded, its evolution semi-group is bounded in any Orlicz space (with bound 
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depending on time). In Section |21 we start looking at solutions of the Boltzmann 
equation. We prove Povzner lemma and several moments estimates in L^, from 
which we deduce the existence and uniqueness part of Theorem 11.21 In Section HI 
we extend the existence result to collision rates depending on the kinetic energy of 
the solution by proving a weak stability result on the basis of (local in time) non- 
concentration estimates obtained by the study of Section |21 to obtain the existence 
part of Theorem 11.41 The uniqueness part of Theorem 11.41 is obtained by proving a 
strong stability result valid for smooth solution. In Section El we study the cooling 
process and prove the remaining parts of Theorem 11.21 and Theorem 11.41 

2 Estimates in Orlicz spaces 

In this section we gather some new functional estimates on the collision operator 
in Orlicz spaces, that will be used in the sequel to obtain (local in time) non- 
concentration estimates. Let us introduce the following decomposition 6 = 6* + 6^ 
of the cross-section b for e G (0, 1): 



where l{-i+£<u.z<i-e} denotes the usual indicator function of the set { — 1 + e < 
u ■ z < 1 — £:}. When no confusion is possible the subscript e shall be omitted. 

In the sequel, A denotes a function strictly increasing, convex satisfying the 
assumptions ()A.1|) . ()A.2|) and ()A.3|) (see the apppendix). This function defines the 
Orlicz space L^(M^), which is a Banach space (see the definition in the appendix). 

2.1 Convolution-like estimates on the gain term 

In this subsection we shall prove convolution-like estimates in Orlicz spaces. These 
estimates extend existing results in Lebesgue spaces: see [201 1211 IHHl 112] in the 
elastic case and JH] in the inelastic case. The proof relies only upon elementary 
tools, essentially Young's inequality, in the spirit of J2]- Moreover it has several 
advantages: its simplicity, the fact that it handles only the dual form of and the 
fact that it is naturally well-suited to deal with Orlicz spaces, since it is based on 
Young's inequality. 

As shown by the formula for the differential of the Orlicz norm in the appendix, 
the crucial quantity to estimate is 



Most of the difficulty is related to the fact that the bilinear operator is not 
bounded because of the term |f — in the collision rate. Nevertheless it is possible 
to prove a compactness-like estimate with respect to this algebraic weight. When 



(2.1) 



h\{£, u; dz) 
bl{S, u; dz) 



b{S, u; dz) l{_i+£<^.^<i_£} 
b{S, u; dz) — bl{S, u; dz) 
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combined with the damping effect of the loss term this estimate shall show that the 
evolution semi-group of the global collision operator is bounded in any Orlicz space. 
Let us state the result 

Theorem 2.1 Assume that B satisfies 111. 5]) - n~Tl\} - lll. l'^) - M. 1 '^] . For any function 
/ G n L^, for any e G (0, 1), there is an explicit constant C^{e) such that 



(2.2) +(2 + 2"«)j£(e)||/||^ / s^'^^\\^,\dv 



LA 



LA 

Remark 2.2 Let us comment on the conclusions of this theorem. 

1. We establish estimates for the quadratic Boltzmann collision operator but 
similar bilinear estimates could be proved under additional assumption on b, namely 
that either no frontal collision occurs, i.e., b{S, u; dz) should vanish for u close to z, 
or no grazing collision occurs, i.e., b{£, ; dz) should vanish for u close to —z. For 
more details on these bilinear estimates and the corresponding assumptions, we refer 
to VJ'Jl where they are proved in Lebesgue spaces in the elastic framework. 

2. Let us emphasize that for 2; ~ (close to sticky collisions), the jacobian of the 
pre-postcollisional change of variable (f,f*) {y',v[) (both velocities at the same 
time) is blowing up. However in our method, we only use the changes of variable 
V v' and v,, — > v' , keeping the other velocity unchanged, and the jacobians of these 
changes of variable remain uniformly bounded as z — * 0. This explains why our 
bounds includes the sticky particules model, and are uniform as z ^ 0. 

3. When A(t) = t^/p, estimate just writes 



(2.3) / Q''ifJ)f'-'dv < Ct{e) ll/ILj \\f\\%+Ue) II/IIlJ 11/1^1^/111 



LP1 



for any e G (0,1) and for some explicit constants C'^{e), jei^) £ (0, 00) with 
jsi^) ~^ when e ^ 0. Although the quantities involved in these previously men- 
tioned papers are slightly different, one can see that estimate \2. (or the L^ version 
of Theorem \2.(^] generalizes \12i Proposition 2.5] to the inelastic Boltzmann opera- 
tor and that it improves 'W, Lemma 4-i] because of the better control of the norm 
\\f\v\ 



l^/^l|L.. 



Let us start with an elementary geometrical lemma that we shall need several 
times in the sequel, in order to justify the change of variables f * v' (keeping v, z 
fixed) and v ^ v' (keeping v^,, z fixed). This lemma is close to the spirit of the proof 
of these changes of variables in the proof of the so-called "cancellation lemma" in 
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Lemma 2.3 For any z E D and 7 G (—1, 1) we define the map 

(2.4) u^w = <^,{u):=u+\u\z, 

its Jacobian function Jz := det {D $2) and the cone := {u G M^\{0}, -u ■ i > 7}. 
Then $^ is a C°°-diff'eomorphism from fi^ onto Qs with 



S = 



7 + 


2; 




(l + 27|z| 


+ 


\z\ 


|2)l/2 



and there exists Cj G (0, 00) such that 

(2.5) C-i < < on fi, 

uniformly with respect to the parameter z E D. 



Proof of Lemma \2.3l We may assume z ^ since otherwise the conclusion is clear. 
Let start proving that $2 is one-to-one on = R^\{R_z). For any x G we 
introduce the decomposition x = Xi z + X2 ■= {xi,X2) such that Xi G M, X2 G M^, 
X2 ■ z = 0. The expression ()2.4|1 then writes equivalently 

Wi = Ui + {ul + \u2\'^y^'^ \z\, W2 = U2. 

For any u,u' G the relation ^^(u) = ^^(m') =: w implies immediately U2 = 
u'2 = W2 and we conclude observing that for any z E D and W2 G the map 

V^u,2,kl : Ml u^i := t^i + {ul + \W2\'^Y^^ \z\ 

is strictly increasing from M onto M if |2;| < 1, from M onto IR+ if |2;| = 1 and W2 0, 
and from IR+ onto IR+ if = 1 and W2 = 0. That proves that $ is one-to-one. 
Moreover, any point u = (^1,^2) G such that -Ui = 7 is mapped to the point 

w = {'-f + \z\,U2) with square norm |wp = 1 + 27|2;| + We conclude that 

$(r2^) = Qs thanks to the homogeneity property ^^(rM) = r^^iiJ-) for any r > 
and u G M^. We next compute D$z(m) = Id + it ^ z and thus J2(m) = 1 + it ■ z 
from which ()2.5j) easily follows. Finally, the fact that $2 is a C°°-diffeomorphism is 
a direct consequence of the local inversion Theorem. □ 

Proof of Theorem \2.1l Let us denote 

ilk)- 

Using the decomposition b = + , we control separately the two terms /* and P 
in the decomposition 

[ Q+{f,f)vif)dv = [ ffMf')\u\b\£,u-dz)dvdv, 

+ / ff,^{f')\u\V{£,u-dz)dvdv,=:l' + r. 
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Using the bound 

\u\ = \v - v^,\ < \v\ + 

we have 

I' < [ {f\v\)fMf')b\S,u;dz)dvdv, 

f{f*\v*\)ip{f') \u\b\S,u;dz)dvdv^ =: !{ + l\. 

'xR^xD 

For the term !{, by applying the Young's inequahty ()A.4|) 

fMf) = II/IIm f ^) vif) < II/IIla a (jjfA + II/IUa a*(^(/')), 



we get 



II < II/IIm / f\v\A(j^] b\£,u-dz)dvdv, 

'xRi^xD VII/IUa, 



+ II/IIla / /|t;|A*(^(/')) b\S, u- dz) dvdv, =: I^, + 

jR^xR^xD 

On the one hand, using 

VxgM+, A(x)<xA'(x), 
which is a trivial consequence of the fact that A(0) = and A' is increasing, we have 

4i<a(^)||/IUl / f^U)dv. 
Jrn 

Holder's inequality in Orlicz spaces ()A.5|) recalled in the appendix then yields 



(2.6) II < a{£) N^* ( A' ( ^ ) ) ||/|Lj ||/|Ua. 



LA. 

On the other hand, using that A*{y) = y (A')"^(y) - A((A')"^(?/)), we get 

11,2 < f f\v\^{f)f'b\£,u-dz)dvdv,. 

We make the change of variables — > v' (while the other integration variables 
are kept fixed) or more precisely \E' : {v,v^,z) {v,i/jy^z{v^), z) with ipv,z{v^) = 
v' = V + 2~^^z{v* — v). Thanks to the truncation ()2.1|) on 6* and Lemma 12.31 
the application \E' is a C°°-diffeomorphism from {{v,v^,,z) G x D, u ■ z ^ 1} 
onto its image and its jacobian Jij, = 2~^ (1 — u ■ z) satisfies |Jv^"'^| < 2^ on 
{(f , f*, z) e M^^ X L), u- z<l-e}. We then get 

lU < f f\v\fv{f) J^H\S,v ~^-l{v'y,dz)dvdv' 

Jr^xR'^xD 
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As previously, Holder's inequality ()A.5|) then yields 

(2.7) ll,<a{S)2^e-'\\fh.N^' {m^)) 

Next, the term is exactly similar to /*, except that one has to use the change 
of variable u — >• f ' = f * + ^z{v — v*) instead of — > v'. Therefore, gathering 
fl2.(ij) . ()2.7j) and the same estimate for I2, we obtain 



N^' i A' ' 



(2.8) /*< 2 a(^) (1 + 2^5-1 

Finally, for the term we can split it as 

r < f f^ip{f')l{u.z<o}\u\b''{S,u;dz) dvdv^ 

+ f f* <^{f') l{«.^>o} \u\ V{£, u; dz) dv dv^ =: /[ + I^. 

For /[, we use Young's inequality ()A.4|) on a; = and y = ^p{f') to obtain 
/[ < / / </)(/*) l{u.z<o} \u\ u] dz) dv dv^ 

JR^ xR^ xD 

ff ip{f') l{u.z<o} \u\ u; dz) dv dv^. 



'xR^xD 

In the second integral we make again the change of variable defined by \E' for which 
there holds < 2^ on the domain of integration because of the truncation 

ii ■ z < 0. We also observe thanks to a direct computation starting from ()1.3|) that 
under the truncation u ■ z < there holds 

\v - v,\ < 2\v' -v\< 2(1 + |f'|)(l + \v\). 

Hence we obtain 

/[ < (l + 2^+^)fsup / b^{S,u;dz))\\fhi [ f ^if) + dv 
< (l + 2^+i)a(f)j,(e)||/|U; / f^if)il + \v\)dv. 



The term I2 is treated similarly using Young's inequality (this time on x = f and 
y = (p{f')) and the change of variable v ^ v' instead of — >• v' . It satisfies therefore 
the same estimate. Thus we obtain the estimate 

(2.9) r< (2 + 2^+2)^(^)^.^(^)11^11 j /^(/)(i + |^|)rf^. 
Defining 

(2.10) Ct{e) = 2(1 + 2^£-i) + (2 + 2''+^)je{e), 

we conclude the proof gathering (j2.8|) and (j2.9j) . □ 
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2.2 Minoration of the loss term 

In this subsection we recall a well-known result about the minoration of the loss 
term Q~ . Let us recall first the following classical estimate. 

Lemma 2.4 For any non-negative measurable function f such that 
(2.11) / fdv = l, [ fvdv = 0, 



we have 

VfGM^, / f,\v -v,\dv,>\v\. 



Proof of Lemma 12. 4[ Use Jensen's inequality 



with the probability measure dfi^: = df*, the measurable function f = w — f* 

and the convex function (p{s) = \s\. □ 

Then the proof of the following proposition is straightforward: 

Proposition 2.5 Assume that B satisfies M.5\) . For a non-negative function f 
satisfying \2.11\) . we have 



(2.12) / Q-if,f)A'(-^) dv>aiS) [ f A' ( 



/ 



LA 



If I dv. 



2.3 Estimate on the global collision operator and a priori 
estimate on the solutions 

Combining Theorem 12 . 1 1 and Proposition 12.51 we get 

Theorem 2.6 Assume that B satisfies M .^) - n~Tl\) - M . 1"^) - M . l'J\i . Let us consider 
a non-negative function f satisfying \2.11]) . Then there is an explicit constant Cs 
depending on the collision rate through the functions a and js such that 



LlWJ IIM. 



More precisely, Cg = a{£) C^{£o), with Eq such that j£{eo) < (2 + 2^+^) ^ 
and where is defined in i2.1(J\} . 



-1 

L} 



Proof of Theorem \2. 61 One just has to combine ()2.2|1 and ()2.12|1 and pick a Eq small 
enough such that 

(2 + 2^+2) \\fhue{eo)<l. 



□ 
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Corollary 2.7 Assume that B satisfies /ll.^) - n~n\) - /ll. r^) - /ll. and let us con- 
sider a solution f G C{[0,T]] L^) to the Boltzmann equation / li.i)) -/ f7T^) associated 
to an initial datum /in G L\ and to the collision rate B. Assume moreover that 
holds and there exists a compact set K C (0, +00) such that 

VtG[0,r], S{t)eK. 

Then, there exists a C'^ , strictly increasing and convex function A satisfying the 
assumptions \A.l}) . \A.^) and \A.3\) (which only depends on fin) and a constant Ct 
(which depends on K, T and B) such that 

snp\\f{t,.)\\LA<CT. 

[0,T] 

Remark 2.8 Let us emphasize that these non-concentration bounds are valid for the 
sticky particules model (in this case they provide an exponentially growing bound in 
for all times). As a particular case we deduce some explicit bounds on the entropy 
when it is finite initially. Moreover, since our bounds are uniform as b ^ ^z=0; '"^e 
also deduce a proof of the sticky particules limit (for a cross-section being a diffuse 
measure converging to a Dirac mass at z = 0) by the Dunford- Pettis Lemma. This 
shows moreover that this limit is not singular. 



Proof of Corollary \2. 71 Since /in G L^(M^), as recalled in the appendix, a refined 
version of the De la Vallee-Poussin Theorem [25, Proposition 1.1.1] (see also [2311211) 
guarantees that there exists a function A satisfying the properties listed in the 
statement of Corollary 12.71 and such that 



/ A(|/in|)rft; < +00. 



Then the norm of / satisfies 



d 
di 



Wfi 



til LA 



A' 



I/I 



LA 



-1 



Q(/,/)A' 



I/I 



LA 



dv 



thanks to Theorem IA.21 and thus using Theorem 12.61 we get 

VtG[0,T], ±\\f^\\^,<Ceit)\\ft\\Ll\\ft\\L-. 

Thanks to the assumptions on B, the constant C£(^t) provided by Theorem 
uniform when the kinetic energy belongs to a compact set. Thus we deduce 



(2.13) 



d 



VtG[0,T], -\\f,\\^,<CK\\ft\\Ll\\ft\\L-- 



IS 



for some explicit constant Ck > depending on K and the collision rate. We 
conclude thanks to a Gronwall argument. □ 
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3 Proof of the Cauchy theorem for non-coupled 
coUision rate 



In this section we fix T^, > and we assume that the colhsion rate B satisfies 

(3.1) B = B(t, u; dz) = \u\ •yit) bit, u; dz), 

where 6 is a probabihty measure on D for any t G [0, T*] and u G M.^ satisfying 

(3.2) Vt G [0,T,], Vm G M^, b{t,u; dz) = b{t, -u; -dz) 
and where 7 satisfies 

(3.3) 0<7(t)<7, on (0,T,). 
3.1 Propagation of moments 

In this subsection we establish several moments estimates which are well known 
for the Boltzmann equation with elastic collision, see [Ul |^ |2S1 and the references 
therein, as well as the recent works jTSl E] for the inelastic case. Let us emphasize 
that these moment estimates are uniform with respect to the normal restitution 
coefficient e or more generally to the support of b{t,u; ■) in D. 

First we give a result of propagation of moments valid for general collision rates 
using a rough version of the Povzner inequality. 

Proposition 3.1 Assume that B satisfies For any < G L^iR'^) 

with q > 2 and T > 0, there exists Ct such that any solution f to the inelastic 
Boltzmann equation M.1]) - I[T^) on [0, T] satisfies, at least formally, 

sup||/(t,-)IUj<CT. 



[0,T] 



Proof of Proposition I3.il We write the proof for the third moment, the general 
moment estimate being similar. For any function : —>■ R_|_ such that \l/(f) := 
^/'(If P) for some function ip : ]R_|_ ]R_|_, the evolution of the associated moment is 
given by 

— / f^idv= I ff^Kq,dvdv^, 

where 

K^ := \ ! (^' + - ^ - ^*) B{t, u- dz). 
^ Jd 

For ip{z) = z^, s > 1, the function ip is super-additive, that is ip{x)+ipiy) < i'ix+y), 
and it is an increasing function. As a consequence, 

^'+^:-^-^, < ^{\v'\^)+^{\v'f)-^{\vf + \vX) 

+H\v\'' + \v.\'') - ^/j{\v\^) - ^i\v,\^) 

< ^(|^|2+|^^|2)_^(|^|2)_^(|^^|2)^ 
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which imphes 

Making the choice ip{x) = x^/^ and using the inequahty 

(a;i/2 + ^i/2)[(^ + ^)3/2_^3/2_^3/2] < C {x'/' + y'/') (x'^'y + xy'^') 

(3.4) < C {2xy + x'/Y^' + x'/y/^) 
for any x,y > 0, we get 

(3.5) ^[ f\vfdv<C-f{t) [ f M\vf\v,f + \v\\v^f)dvdv,, 

and we conclude thanks to a Gronwall argument. □ 

Finally we give a much more precise result on the evolution of moments in the 
case when assumption H4 is made. On the one hand, we state uniform in time 
propagation of algebraic moments (as introduced in jHU 12 CHI) and exponential 
moments (for which the first results were obtained in [Hj). On the other hand, 
we prove appearance of some exponential moments (while appearance of algebraic 
moments was initiated in EHl EH]) using carefully estimates developed in jH]. 
These estimates may be seen as a priori bounds, but in fact, by the bootstrap 
argument introduced in they can be obtained a posteriori for any solution 
given by the existence part of Theorem 11.21 and Theorem 11.41 

Proposition 3.2 We make the assumption H4 on B. A solution f to the inelastic 
Boltzman equation / ti.i)) -/ f7T^) on [0, Tc) satisfies the additional moment properties: 

(i) For any s > 2, there exists Cs > such that 

(3.6) sup ||/(t,.)lUi <max{||/i,|Ui,Cj. 

te[o,Tc) 

(ii) If fine'' ^"^^ e Li(R^) /orr > andrj G (0,2], there exists Ci,r' > 0, such that 

(3.7) sup / f{t,v)e'''\''\" dv <Ci. 

(Hi) For any rj G (0, 1/2) and t G (0, T^) there exists a^, G (0, oo) such that 

(3.8) sup [ /(t,t;)e""l^l''rff < C^. 

Let us emphasize that none of these constants depends on the inelasticity coefficient 
e (so that the estimates are uniform with respect to the inelasticity of the Boltzmann 
operator) and that the constant Cs,ar,, may depend on /in only through its kinetic 
energy Ei^- 
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Remark 3.3 The proof of (i) is very classical for the elastic Boltzmann equation 
(73) / and it has been extended to the inelastic operator in JW^ . Estimate (ii) 
has been proved in '^61 for the elastic Boltzmann equation and it has been generalized 
in ^ to the (stationary) inelastic Boltzmann equation. We refer to \31\. 1^ 
for development around the Povzner inequalities. Since (ii) is a straightforward 
consequence of the Povzner inequality proved in we just have to prove (Hi). 
Nevertheless, since the proof of (Hi) requires some tools and notations introduced 
in we begin (step 1 and step 2) by briefly presenting the proof of (ii). Let 

us emphasize again that (Hi) is new even for the elastic equation. In the elastic 
framework, an extension of (Hi) to hard potentials with cutoff has been used recently 
in the proof of the exponential return to equilibrium with explicit rate for initial data 
with finite mass and energy, see U^ . 

Proof of Proposition VA.'A Let us define 

nip := / / |f dv. 



Step 1. Differential inequalities on tlie moments. Taking ip(x) = x^l"^ and B of tlie 
above form, there holds 



(3.9) —mp = / Q[f , f)\vf^ dv = a{E) \ f f^\v - Kp{v,v^) dv dv^ 
where 



(3.10) Kp{v,V,):=l I (l^ff 



2 7g]v-i * a{S) 

From [HI Lemma 1, Corollary 3], there holds 
(3.11) Kp{v,v,)<-ip{\v\^+\v,\') 
where (7p)p=3/2,2,... is a decreasing sequence of real numbers such that 

4 



U, |2^P_ |^|2P_ 1^ |2P 



(3.12) 0<7p<minn 



p+ 1 



(notice that the assumptions |Hl (2.11)-(2.12)-(2.13)] are satisfied under our assump- 
tions on the collision rate). Let us emphasize that the estimate ()3.11|) does not de- 
pend on the inelasticity coefficient e{£, \u\). Then, from [3 Lemma 2 and Lemma 3], 
we have 

(3.13) ^/ QUJ)\v?''dv<^pSp-{l-^p)mp+y^ 



a{8) 



with 



k=i ^ ^ 



("^fc+i/2 ^p-k + 'Tifc mp_fc+i/2), 
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where kp := [{p + l)/2] is the integer part of (j9 + l)/2 and stands for the 

binomial coefficient. Gathering (j3.9|) and (j3.13p . we get 

(3.14) ^rup < a{£) {-fp Sp - {1 - -fp) mp+1/2) Vp = 3/2,2,... 

By Holder's inequality and the conservation of mass. 



nip < mp+1/2 



and, by [SI Lemma 4], for any a > 1, there exists A > such that 

Sp < AT{ap + a/2 + l)Zp 

with ^ 

k^'^X^'''^"^ ^f-'^+Vs}' ■■= Y{ap + 1/2)- 
We may then rewrite ()3.14|1 as 

,3.15, f (a - (1 - r(.p + 

for any p = 3/2, 2, . . . On the one hand, from (|3.12p . there exists A' such that 



On the other hand, thanks to Stirling's formula n\ e " \/2wri when n ^ 00 

and the estimate ()3.12|) . there exists A' > such that 

(3.17) (l-7p)r(ap+l/2)i/2p> Vp = 3/2,2,... 
Gathering ()3.15p . ()3.16|) and ()3.17|) . we obtain the differential inequality 

(3.18) ^ < a{S) Zp - A"p'^'^ zl+^'^P) 
for any p = 3/2, 2, . . . 

Step 2. Proof of l\3.7]) . On the one hand, we remark, by an induction argument, 
that taking po := max{3/2, {2A'/A"y}, the sequence of functions Zp := is a 
sequence of supersolutions of ()3.18|) for any x > and for p > Pq. On the other 
hand, choosing xq large enough, which may depend on po, with have from (i) that 
the sequence of functions := is a sequence of supersolutions of ()3.18|) for any 
X > Xq and for p G {3/2, . . . ,po}. As a consequence, since Zp for p = 0, 1/2, 1 are 
bounded by ||/in||L^! we have proved that there exists xq such that the set 

(3.19) := |z = (zp); ZpKx^ Wpe^N 
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is invariant under the flow generated by the Boltzmann equation for any x > Xq: if 
f{ti) e then 7(^2) e Ca; for any ^2 > ^i- We set a := 2/r] > 1. Noticing that 



(3.20) / /(t;)e^l'^l''rfi; = V^m,,/2 

we get, from the assumption made on /in, that 

A;' 

mfe/a(0)<Co^ VA;eN. 

Since we may assume r G (0, 1], the function y Cq r(?/ + l)r~^ is increasing, and 
we deduce by Holder's inequality that for any p 

^M<CJ^^<CJ-^^^ with £,:=M + 1. 
From the definition of Zp we deduce 

(3.21) ,^(0)<C,^?ii^|±i)<.; 

for any p and for some constant Xi G (0,oo). Choosing x := max{xo,a;i} we get 
from and ()H.21|1 that Zp{t) < x^ E [0,Tc) for any p. Therefore, we have 

mp{t) <V{ap+ 1/2) x^ = 3/2, 2, . . . , VtG[0,Tc). 

The function y ^ r(?/ + l/2) x^ being increasing, we deduce from Holder's inequality 
that for any k eW that mk/ait) < V{ap + 1/2) x^ < V{k + a/2 + 1/2) x^/'^+^Z^ ^j^j^ 
p := [2A;/a]/2 + 1/2. For r' < 2x^^/''(l + a)~^ we conclude 

Vt G [0, TJ / /(t, i;) e'-' 1^1" di; < V ^(fc + a/2 + 1/2) ^./.^^/^ 
A- ^ A;! 

from which ()3.7j) follows. 

Step 3. Proof of (13.8)) . Let us fix r G (0,Tc). We claim that there exists x large 
enough and some increasing sequence of times itp)p>pg which are bounded by r such 
that for any p 

(3.22) VtG[tp,Tj Zp{t)<xP. 

We already know by classical arguments (see [HH EHl) that for pq (defined at the 
beginning of Step 2) there exists Xi, larger than xq defined in ()3.19p . such that ()3.22|) 
holds for any p < Po and tp = r/2. We then argue by induction, assuming that for 
p > Po there holds: 

(3.23) Zk<x'' on [tp_i/2,Tc) \/k<p-l/2 

(3.24) Zp>xP on [tp_i/2,tp), 
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for some a; > Xi to be defined. If (I3.24j) does not hold, there is nothing to prove 
thanks to Step 2. Gathering (j3.23p . (j3.24|) with (j3.18|) we get from the definition of 
Po and the fact that S{t) G [S{t),S{0)] so that (y{S) > ao > 

(3.25) ^^-"°T^"^''^^'^'' (Vi/2,g- 

Integrating this differential inequality we obtain 

Defining [tp) in the following way: 



r TV I sr^ 1 0/9 

to:=-, tp:=tp_^i2^-- , s„ :=}_^pi-'^/2 



p=0 



and defining X2 := (8 s^)^/ (A" ao t)^ we have then proved Zp{tp) < and therefore 
Zp{t) < for any t > (tp,Tc) with x = max{xi,X2} thanks to Step 2. Setting 
a := 2/7] > A [r] < 1/2) we have 

00 

(3.26) J2^i+k/2<r 

k=0 

and we conclude as in the end of Step 2. □ 

3.2 Stability estimate in L\ and proof of the uniqueness part 
of Theorem 11.21 

Proposition 3.4 Assume that B satisfies For any two solutions f and 

g of the inelastic Boltzmann equation / li.i)) -/ rO)) on [0,T] (T < T^) we have 

(3.27) ^/ \f-g\{l + \v\^)dv<C^J {f+g){l + \v\')dv j \f -g\{l + \v\^) dv. 



We deduce that there is Ct > depending on B and sup^gjg.T] 11/ + fi'lU;^ such that 
VtG [0,T], \\ft-gt\\Ll < ll/in-^7in||L^e^-*. 

In particular, there exists at most one solution to the Cauchy problem for the inelastic 
Boltzmann equation in C{}fi,T]] L\) nL-'^(0,T; L\). 

Proof of Proposition We multiply the equation satisfied by (/ — g) by (pit, y) = 
sgn(/(t, y) — g(t, y)) k , where k = {1 + |f P). Using the chain rule ()1.3H) . we get for 
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alH > 

d f . „ , . , 1 



dt 



\f-g\kdv = I [if - g)9* + fif* - g*)] 

{(j)' + (j)'^ — (j) — (j)^) B{t, u] dz) dv^ dv 

(0' + 0* — — 0*) u\ dz) dv^ dv 
< I [ \f-9\{f* + 9.) 

^ Jr^xR^^xD 

{k' + k'^ — k + k^) B{t, u] dz) dv^ dv, 



where we have just use the symmetry hypothesis (jS.ljl . ()3.2|1 on B and a change of 
variable — Then, thanks to the bounds (j3.ip . (j3.3p we deduce 

-j: I \f-9\kdv < 7* /" |m| 1/ - 5f| (/* + 5f*) /i;*(it;*(if 
< 7* / \f-9\kdv I {f* + 9*)kl'^ dv, 

which yields the differential inequality ()3.27|) . The end of the proof is straightforward 
by a Gronwall argument. □ 

The uniqueness in C([0,T); L\) fl L^(0,T; L3) as stated in Theorem 11.21 is given 
by Proposition 13.41 

3.3 Sketch of the proof of the existence part of Theorem 11.21 

As for the existence part, we briefly sketch the proof. We follow a method introduced 
in |31j and developed in ^Tj. We split the proof into three steps. 

Step 1. Let us first consider an initial datum /in satisfying ()1.10|) with g = 4 and 
let us define the truncated collision rates Bn = -Bl|„|<n. The associated collision 
operators Qn are bounded in any L^, q > 1, and are Lipschitz in on any bounded 
subset of L^. Therefore following a classical argument from Arkeryd, see j2], we 
can use the Banach fixed point Theorem and obtain the existence of a solution 
< /„ e C{[0,T];Ll) n L°°(0,r;Li) for any T > 0, to the associated Boltzmann 
equation (IIII1)-(IE21), which satisfies (fO^ - (fn?H|l . 

Step 2. From Proposition 13. 11 for any T > 0, there exists Ct such that 

sup ll/nlUi < Ct. 

[0,T] 

Moreover, coming back to the proof of Proposition 13.41 (see also the first step in the 
proof of [T71 Theorem 2.6]), we may establish the differential inequality 

d C2 

^ll/n ~ /mllii < Ci Wfn + frnWhl Wfn — /mlU^ + ~ Wfn + /mHii 
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for any integers m > n. Gathering these two informations we easily deduce that (/„) 
is a Cauchy sequence in C([0, T]; L^) for any T > 0. Denoting by / G C([0, T]; Ll) fl 
L^{0,T] L\) its hmit, we obtain that / is a solution to the Boltzmann equation 
()l.H) - p.2|) associated to the collision rate B and the initial datum /in by passing to 
the limit in the weak formulation ()1.30|) of the Boltzmann equation written for /„. 

Step 3. When the initial datum satisfies (jl.iup with g = 3 we introduce the 
sequence of initial data /m/ := fin'^{\v\<e}- Since /in/ G L\, the preceding step give 
the existence of a sequence of solutions fe G C{[0,T]; L^) fl L°°{0,T; Ll) for any 
T > to the Boltzmann equation ()1.H) - ()1.2|) associated to the initial datum /in/. 
From Proposition 13. H for any T > 0, there exists Ct such that 

sup ||/f||L.i < Ct- 

[0,T] 

Thanks to (j3.27|) we establish that (/^) is a Cauchy sequence in C([0, T]; Ll) and we 
conclude as before. □ 

Remark 3.5 Note here that an alternative path to the proof of existence could have 
been the use of the result of propagation of Orlicz norm which shows that the solution 
is uniformly bounded for t E [0,T] in a certain Orlicz space. Together with the 
propagation of moments and Dunford- Pettis Lemma, it would yield the existence of a 
solution by classical approximation arguments and weak stability results as presented 
below. More generally the propagation of Orlicz norm by the collision operator can 
be seen as a new tool (as well as a clarification) for the theory of solutions to the 
spatially homogeneous Boltzmann equation with no entropy bound, as in the inelastic 
case, or in the elastic case when the initial datum has infinite entropy, see also f^Ul^ 
where other strategies of proof are presented. 

4 Proof of the Cauchy theorem for coupled coUi- 
sion rate 

4.1 Weak stability and proof of the existence part of Theo- 
rem 11.41 

Proposition 4.1 Consider a sequence Bn = Bn{t, u; dz) of collision rates satisfying 
the structure conditions \3.1\) - KS~^) and the uniform bound 

0<7n(t)<7T VtG[0,T], VnGN*, 

and let us denote by fn G C([0,T); Ll) fl L°°(0,T; Ll) the solution associated to Bn 
thanks to the existence result of the preceding section ( existence and uniqueness part 
of Theorem \l.^ and Remark \l.'^ Ath point). Assume furthermore that {fn) belongs 
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to a weak compact set of L^{{0,T) x M^) and that there exists a collision rate B 
satisfying / rO)) -/ f3^) -/ r73|) and such that for any ijj e Cc(K^) 

7n — ^ 7 and / 'ip{v') bn{t,u; dz) — > / 'ip{v')b{t,u;dz) a.e. 

J D J D 

Then there exists a function f G C([0, T); L^) fl L°°(0, T; L3) anc? a subsequence fn,. 
such that 

fn, ^ / weakly m ^^((O, T) x M^), 
and f is a solution to the Boltzmann equation il.l}) - l[T^) associated to B. 

Such a stability/compactness result is very classical and we refer to [3 E] for 
its proof. 

Proof of the existence part of Theorem \1.4l We assume without restriction that 
there exists a decreasing function ao such that a < ao on [0,£^in]. We proceed in 
three steps. 

Step 1. We start with some a priorihounds. We set := H/Hi^i. From the Povzner 
inequahty (|3.5|) (with 7(t) = a{£{t))) and the dissipation of energy equation (|1.8|) . 
we have 

(4.1) jYs < Ci ao{S) Y,, Y,{0)Y,{f,^) 
and 

(4.2) j^S > -Ci ao{S) Ys, S{0) = S,^, 

for some constant Ci (which depends on Si^). There exists such that any solution 
{Y3,£) to the above differential inequalities system is defined on [0,T*] and satisfies 

(4.3) sup Y^it) < 2F3(/in), inf S{t) > 

[0,T.] [0,T.] 

More precisely, we choose T^, such that 

Ci ao(^in/2) < r3(/in) and Ci «o(^in/2)2 F3(/in) < 

in such a way that if (F3, satisfies V3 < 2 Y-^^fi^ and ()4.2|) on (0, T*) or if (Kj, £") 
satisfies E > E;„/2 and dHH) on (0, T^) then (jOl) holds. We introduce 

X:=\Ee C([0,T,]), E-J2 < E{t) < E,^ on (0,T,)}. 

Step 2. Let us consider a function £^1 G X and define B2(t,u; dz) := B{Ei(t),u; dz). 
From assumptions p.5|) - (P~TT|) - ()1.12'|) - ()1.13l we may write 

B2{t, u; dz) = \u\ 72(t) b2(t, u; dz) 

where 62 is a probability measure and 72 (t) satisfies 

72 (t) = a(^i(t)) < «o(W2) < +00 Vt G [0,T,]. 
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Thanks to Theorem ll.2l there exists a unique solution /2 G C([0, T*]; L2)nL°°(0, T^,; Lg) 
to the Boltzmann equation (jl.l|) - (jl.2|) associated to the colhsion rate B2 and we set 
£2 '■= £{.12)- In such a way we have defined a map $ : X ^ X, $(£^1) = £^2- 

In order to apply the Schauder fixed point Theorem, we aim to prove that $ 
is continuous and compact from X to X. Consider a sequence of X which 
uniformly converges to £1. Since belongs to the compact set [£^in/2, i^jj for any 
n and any t G [0,T*], we deduce by applying CoroUarv 12.71 to the sequence (/^) 
associated to Bl^it^ u; dz)B{£^{t),u; dz) that 

(4.4) Vn>0, sup / A(/2"(t,t;))rft; < C2, 

for a superlinear function A and a constant C2 > 0. Moreover, from Proposition 13. II 
we have 

(4.5) Vn>0, sup / f^{t,v)\v\^dv <C3 

for some constant C3 > 0. 

On the one hand, gathering ()4.4|) . ()4.5|) and using the Dunford- Pettis Lemma, we 
obtain that (/^) belongs to a weak compact set of L^((0,T*) x M^). Propositon 14.11 
then implies that there exists /2 G C{[0,T^]; L^) fl L°°(0, T^,; L3) such that, up to a 
subsequence, ^ /2 weakly in L^(0,T; L^) and /2 is a solution to the Boltzmann 
equation associated to B2{t, u; dz) = B{£i(t), u; dz). Since this limit is unique by the 
previous study, the whole sequence (/2 ) converges weakly to /2, and in particular 

(4.6) 8^ £2 weakly in L^(0,T) 

where £2 is the kinetic energy of f2- 
On the other hand, there holds 

^£2 = - [ f2 fl \n? A(f r, u) dvdv. =: -D^. 

Since A(^f,u) < a(^i")/4 < ao(^in/2)/4, we deduce from (jHUI) that is bounded 
in L°°(0,T) which in turn implies 

(4.7) ||i^2 IIh/i.°°(o,t) 

From Ascoli's Theorem we infer that the sequence {£2) belongs to a compact set 
of C([0,T]). Since the cluster points for the uniform norm are included in the set 
of cluster points for the L^(0,T) weak topology, it then follows from (j4.6|) that 
$(£") = £{f2) converges to ^{£1) = £{f2) for the uniform norm on C([0, T]), which 
ends the proof of the continuity of Of course, the a priori bound ()4.7|) and Ascoli's 
Theorem also imply that $ is a compact map on X. We may thus use the Schauder 
fixed point Theorem to conclude to the existence of at least one £ & X such that 
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$(^) = S. Then, the solution / G C([0, T,]; L^) n L~(0,T,;Li) to the Boltzmann 
equation associated to B{t,u] dz) := B{E{t),u; dz) satisfies 

/ f{t,v)\v\^dv = ^{B){t) = B{t) 

Jrn 

and therefore / is a solution to the Boltzmann equation associated to B in 
C([0,T,];Li)nL-(0,T,;Lj). 

Step 3. We then consider the class of solution / : (0,Ti) L\ such that / G 
C([0,T]; L\)P\L°°{Q,T] L3) for any T G (0,Ti), £ is decreasing, / is mass conserving. 
By Zorn's Lemma, there exists a maximal interval [0, Tc) such that 

(Tc < 00 and when t Tc) or Tc = +cxd. 

In order to end the proof, the only thing one has to remark is that if Tc < +00 
and lim £{t) =£c>0, then lim Y^it) < 00 (by dHH)) so that / G C{[0,Tc];Ll) n 

L°°{0,Tc; L3) and we may extends the solution / to a larger time interval. □ 

4.2 Strong stability and uniqueness part of Theorem 11.41 

In this subsection we give a quantitative stability result in strong sense, under the 
additional assumption of some smoothness on the initial datum and the collision 
rate. Let us first prove a simple result of propagation of the total variation of the 
distribution. 

Proposition 4.2 Let B be a collision rate satisfying assumptions ^^H})-^^^-^^^ 

and < /in G BV^HLl an initial datum. Then there exists Ct,, depending on 7,, and 
ll/inllij, such that any solution f G C{[0,T^], L2) fl L°°{0,T^,, L^) to the Boltzmann 
equation constructed in the previous step satisfies 

VtG [0,T,], WMIbv, < Wfinhne''^''. 

Proof of Proposition \4. 'A The proof is based on the same kind of Povzner inequality 
as above. Let us first prove the estimate by a priori approach, for the sake of 
clearness. We have the following formula for the differential of Q: 

V.Q(/,/) = Q(V./,/)+Q(/,V./). 

This property is proved in the elastic case in [3H1 but it is strictly related to the 
invariance property of the collision operator 

rhQ(/,/)=Q(r,/,r,/) 

where the translation operator is defined by 

Vt; G M^, Tt,f{v) = f{v-h). 
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It is easily seen that it remains true in the inelastic case under our assumptions. 
The propagation of the norm has already been established. Then we estimate 
the time derivative of the L\ norm of the gradient along the flow: 



|l|v./tlLi 



/ (V,/,) (1 + \vf) sgn(V,/)' + (1 + \vX) sgn(V,/): 



-(1 + \vf) sgn(V./) - (1 + Iv^l") sgn(V./). 



B dv dv* 



< 



'xRf^ xD 



[l + If J 



/iv.M (i + ifr) + (i+Kn-(i + ifn 



Bdvdv, + A^,\\ft{l + \v\'')\\L^ ||V,/(1 



< c'||/,|Ui||v./|Ui 

using a Povzner inequahty as in ()3.4p . This shows the a priori propagation of the 
BVi norm by a Gronwall argument. 

Now let us explain how to obtain the same estimate by a posteriori approach. 
First concerning the a posteriori propagation of the L\ norm, it is similar to the 
method in [3T] and does not lead to any difficulty. Concerning the propagation of 
BVi norm, we look at some "discretized derivative". Let us denote k = sgn(r/i/ — 
/) (1 + |f 1^). We can compute by the chain rule the following time derivative (using 
the invariance property of the collision operator) 



d 



j 



'xD 



{ThfThf* - //*) [k' - k]Bdv dv^ 
{Thf - /)/* [fc' + fc* - - k^ B dv dv^ 



1 

+ 2 



{Thf - f){Thf* - /*) [k' + k'^-k- k^] B dv dv^ 



< 



'xR"xD 

\Tkf-f\f. 



'xRJVxD 



B dv dvt: 



+ - 



'xR'^xD 



l^hf - f\\rhf* - f* 



.'1 4 



|4 



B dv dvt,. 



Then using the same rough Povzner inequality as in the proof of Proposition 13.11 
we have 



v'\^ + \v:r + \vr + \v. 



J |4 



Hence we deduce that 



\v — vJ < C 



(l + |i;|^)(l + |t;,|^) + (l + |f.|^)(l + |f|^) 



d 



^Ik/i/t - MIli < c't* hhft - MIli II/IIli + llThft - MIli 
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and for \h\ < 1, we deduce 

^Ikh/t - MIli <Ci* Wuft - /tllLill/IUi- 

By a Gronwall argument it shows for any < 1 that 

Vt G [0,T,], ||r,/i - fthi < ||r,/i, - /inllLje^-** 

for a constant Ct, depending on 7* and sup^^jQ ll/tlUi- By dividing by h and 
letting h goes to 0, we conclude that 

Wt e [0,T,], WVJtWMl < llV^/inllMie^^** 

which ends the proof. □ 

Assume now that the collision rate satisfies (jl.5|) - (frTT|) - (jl.l2|) - (jl.l3p plus the 
additional assumption HI. Let us take /in G -BV4 fl Ll and let us consider two 
solutions f,gE C{[0,Tc]; L^) fl L°°(0,T; L3) constructed by the previous steps. For 
these two solutions the function e{S) is locally Lipschitz, so is the function A{S) and 
the differential equation (jl.Sp satisfied by S{ft) on [0,T*] implies that it is bounded 
from below on this interval. Thus thanks to the continuity of a, the assumptions 
of Proposition 14.21 are satisfied, and thus the 5V4 norm is bounded on any time 
interval [0,T,] C [0, . 

Proposition 4.3 Let B be a collision rate satisfying I^l.^} - ^~T1\) - I^l.l'^) - I^l.l^} plus 
the additionnal assumption HI. Let f,g G C([0, T^,]; L2) fl L°°{0,T^,; Ll) be two 
solutions with mass 1 and momentum 0, with initial data fin and g-m, (^nd such that 
S{f{t, .)), S{g(t, .)) E K on [0, T^] with K a compact of (0, +00) and 

VtG[0,T,], \\f{t,.)\\BvA9{t.-)\\BV,<CT,. 

Then there is a constant C^^ depending on B, K and Ct^., such that 
Vt G [0, T,], Wfit, .) - g{t, < ll/in - ^?in|Ui e^-. 

We need the following geometrical lemma which is a more accurate version 
of Lemma 12.31 when the collision process is of the generalized visco-elastic type 
(I1.14I1.15|1 . 

Lemma 4.4 For any e G (0, 1] and a G S^^^ we define 

(4.8) 0: = 0;,,, : M^, v,^v' = v + ^ ^^{v, - v) 

(4.9) 0e = 0e,..,a : ^ M^, v^v' = v, + ^^<^rAv-v,), r, = l^, 

4 3 — e 

(where was defined in Lemma W^) and the J acobian functions J* = det {D (pl^y^t^), 

Je = det (D0e,^^,<^). 
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Then for any 7 G (—1, 1), (p* defines a C°° -diff'eomorphism from v + onto 
V + Quj'i-r) with Lj*{'y) = ((1 + 7)/2)^/^ and (pe defines a C°° -diffeomorphism from 
+ ri-y onto + ^^(^^(7) with uJe{'j) = (7 + re)/(l + 2'jri, + r^Y^"^ . Moreover, there 
exists G (0, 00) such that 

(4.10) C:;^\v~v,\ < \(P^{v)-v,\ <2\v~v,\, 

(4.11) \<P:\v') - (p;,\v')\ < a, |e' - e| \v' - v,\, 

(4.12) I Jel < C7, \J-^\ < C^, \J;' - Jf\ < \e' - e| 

on v^ + Q.y uniformly with respect to the parameters e, e' G [0, 1], cr G G M^. 

The same estimates hold for 0*. 

Finally, for any e, e' G [0, 1], cr G S^^^, G and t G [0, 1] i/iere /ioWs 

(4.13) t0;i + (i_t)0-i = 0-,i 

/or some e" into the segment with extremal points e and e' . The same result holds 
for (PI. 

Proof of Lemma 14.41 We only establish the result for the function (pe, since the 
proof for (pl is similar (and even simpler). First, ()4.12j) and the fact that (pe defines 
a C°°-diffeomorphism from + onto i\ + ^wei^y) come straightforwardly from 
Lemma 12.31 and its proof. 

Second, for any z E D, \ = \u + \u\ z\ < 2 \u\ and 

l^ziu)]"^ > \u\^ + 2\u\z-u + |up \z\^ > |np (1 - 7^) 

for any u G M^, u ■ z That proves 1)4.101) . 

Third, using the notation of Lemma IT^ we write ^~l{w) = {ip~l^{wi),W2) for 
any w = (^1,^2), Wi G M, W2 G M^, W2 ■ <J = 0. The map {ui,r) ip^^^riui) is 
smooth and has positive partial derivatives on Mx [0, 1] if 7^ and on (0, 00) x [0, 1] 
if W2 = 0. On the one hand, we deduce that (wi,r) {p^l^{wi) is smooth and 
increasing in both variables and that the same holds for 

4 

(wi,e) ^ j^VZirS'^i)- 

The intermediate values Theorem then implies that for any e < e' G [0, 1], t G [0, 1] 
there holds 

4 4 4 

^ JZr^ ^^IreM + (1 - ^) Vll,r,, iWl) = J—^ ip;,l,^„ (Wi) 

for some e" G [e, e'] from which ()4.13|1 follows. 

On the other hand, r h-^ ^^^{w) is smooth for any w G S^~^\{— a} and therefore 
there exists such that |$~j;(w) — $~j:(w) | < \r'—r\ uniformly for any w G E>^~^, 
w ■ a > 'J. Thanks to the homogeneity property $^^(A w) = A ^~^{w) we deduce 

\KaM - = MKli^) - Kli^)\ < c., \r' - r\ \wl 
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from which (j4.1H) follows. □ 
Proof of Proposition 14.31 Let us denote Qf (resp. Qg) the collision operator with 
collision rate associated with S = £{f) (resp. £ = £{g)), D := f — g, S := f + g 
and A; := (1 + {vl"^) sgn(D). The evolution equation on D writes 

= i [QfiD,S) + QfiS,D)] + [Qfig,g) - Qgig,g)] 
and thus the time derivative of the L\ norm of D is 



liPii.. 



SD, 



^«(/)) + ^K,e(/))-^-^* 



+ 
+ 
+ 



gg* 



\u\ 



gg* 



gg* 

xRJ^xgJv-i 
: J1 + J2 + /3 + /4, 



if)) 



\u \ b£(^f ) dv dv^ da 



dv dv* da 



\u \ J J dv dv^ da 



\u \ b^f^j-j dv dv* da 



the subscripts recalling that the post-collisional velocities ^^e(j), Ke(f) K{g) 
fined by ()1.28j) depend on the choice of the normal restitution coefficient e and 
thus on the kinetic energies £{f) and £{g)- Here we have set bsi^x) = b{£,x) and 
b^{x) = b{£, x) l±a;>o and for the sake of brevity we just write e{h) instead of e{£{h)) 
for any function h & L\. 

The first term is easily dealt with by the same arguments as in the non-coupled 
case: 



h < 



'xR^xS^-i 



S \D*\ (1 + \v\'') \u\ beif) dvdv* da < «(£:(/)) \\S\\^ \\f - gh^. 



Using \u\ {\k\ + \k{v'^(^g^)\) < 2 (1 + \v\^f/'^{l + |^;*^^/^ the second term h 



IS 



controlled by 

l2<2 - 6£(g)||Li(5^-i) \\g\\li- 

Using now the locally Lipschitz assumption ()1.16|1 and the fact that £{f),£{g) G K 
we get for some constant Ck depending on b and K: 



h < Ck \£{f) - £{g)\ \\g\\% < Ck \\f - ghi Ml- 

As for the third term J3, we use twice the change of variable v ^ v' 
f*, a fixed and e = e(/) or e = e{g)- We get 



(f ) with 



'xS'v-i Jo. 



g* k' J) dv' dv* da 



(/) 



'xS'v-i Jo 



g* k' G{(t)^fl^^) J^J) dv' dv* da, 



<9) 
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where we have introduced the notations G{w) := |w* — w\ g{w) for any w G and 
Oe = + ftui^(o). Without restriction we may assume e(/) < e{g) and therefore 
Oe{g) C Oe(f) since e i— > u;e(0) is an increasing function. We then spht J3 as 



k' G {(/) ,).) J ,\ btr. dv' dv^ da 



+ 
+ 



g* k' 



e{f)J ^e{f) 

^^^e{g)) '^^ 



-^3,1 + -^3,2 + -^3,3- 



g* k' G'(0^(^^)) - J^J) rft;' dv^ da 



For the first term J31 we use the backward change of variables v' ^ v = (p^^^^^v') 
and we get 



h,i= I I g^k{v[,..)Ghsi^f)lo<u.a<r,dvdv^da 

with 7] := i^~(^f^ o i^e(g)(0). By inspection, the functions (e,7) G [0,1] x (—1,1] ^ 
ujf;{'y),uj~^{'~f) G (—1,1] are smooth with respect to both variables. From this 
smoothness and the fact that cu'^ o uje{0) = we deduce \uj'^^ o ^^£(0)1 < C |e — e'| 
for any e, e' G [0, 1] and for some constant C G (0, 00). As a consequence, thanks to 
the Lipschitz assumption ()1.17|) . we obtain 



da > dv dv^ 



h,l < II^IIl- / g (1 + \v\Y g^ (1 + \v^\y <^ I J-C(e(g)-e(/))<«-<7<0 ' 

Jr^xr^ Lj§jv-i 
< C\\g\\l,\e{f)-e{g)\<C\\g\\lAf-g\\Li. 

For the term 73^2, using the estimate ()4.12|) and the Lipschitz assumption ()1.17p . 
we get 

' <C\eif)-eig)\<CK\\f-g\\Ll■ 



J-^ _ 7-1 



Then doing the backward change of variable v' v = y^^^-^^ 
Je(/) is bounded on {u, u ■ a >0} thanks to (j4.1(J|) . we get 

13,2 < Ck hWli- 



'v') and observing that 



We now aim to prove that for any functions /, g which energies £f and £g be- 
longing to a compact K C (0, 00) there exists a constant Ck such that the following 
functional inequality holds 



(4.14) 



4,3 < CxWf -gWa WaWhi M 



BVa- 
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Let us first assume that / and g are smooth functions, say f,gE From 
(lOTIl and we have 



< 



< 



dt 



< C \e{f) - e{g)\\v' - v\ f 

Jo 



dt 



with et G [e(/), e{g)]. Since then Oe(g) C Oet for any t G [0, 1], we deduce 

•1 



h3 < C \e{f) - e{g)\ 



JM^xs^-i 



5'* l^'l I'w' 



Using finally the backward change of variable v' ^ v 
bound ()4.12j) on the Jacobian Je^ on t>* + Qq we get 



VwG{(f)^^^{v')) dv' dv dadt. 
= and the uniform 



43 < C |e(/) - e{g)\ H^H^i ||^||w4- 

Therefore we obtain ()4.14|) for smooth functions. When f,g& BV^ we argue by 
density, introducing two sequences of smooth functions (/„) and ((7„) which converge 
respectively to / and g in and are bounded in -BV4, we pass to the limit n — > cxd 
in the functionnal inequality ()4.14|1 written for the functions and gn- We then 
easily conclude that ()4.14p also holds for / and g. 

The term 4 can be dealt with similarly to the term I^. Collecting all the esti- 
mates we thus get 

-j^Wft - 9t\\L\ < Wft - gthl 

where C^^ depends on K, b and on some uniform bounds on and ||fi'||_BV4- This 

concludes the proof by a Gronwall argument. □ 

The uniqueness part of Theorem 11.41 follows straightforwardly from Proposi- 
tion |01 and the discussion made just before its statement. 



5 Study of the cooling process 

In this section we prove the cooling asymptotic as stated in point (ii) of Theorem 
11.21 and points (iii), (iv), (v) of Theorem 11.41 We first prove the collapse of the 
distribution function in the sense of weak * convergence to the Dirac mass in the 
set of measures. 

Proposition 5.1 Let G (0, +00] be the time of life of the solution. Under the 
sole additional assumption H2, there holds 

(5.1) f{t,.) 5y=o weakly* m M\R^). 

t — >Tc 
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Proof of Proposition 15. Jl We split the proof in two steps. 

Step 1. Assume first that £ ^ when t — >■ Tc. This is always the case when 
Tc < +00 (since the convergence to of the kinetic energy follows from the existence 
proof in this case) and it will be established under additional assumptions on B when 
Tc = +00 but it probably holds true under the sole assumption H2 in this case as 
well. For any < y9 G P(]R^\{0}), there exists r > such that (f = on D{0, r) and 
then, there exists = C^{r, Hv^Hoo) such that \^{y)\ < C^p |f p. As a consequence, 



I fipdv<Cp£it)^0, 



from which we deduce that any weak * limit yU of / in satisfies supp yU C {0}. 
Therefore, fl5.1|l follows using the conservations ()l.H2j) and the energy bound (jl.HHj) . 

Step 2. Assume next that £ — > £qo > (and thus also Tc = +oo). Then for a 
fixed time T > and for any non-negative sequence (in) increasing and going to 
+00, there exists a subsequence (tru.) and a measure G L^{0,T;M^) such that 
the sequence fkit,v) := f(tn^ +t,v) satisfies 

(5.2) fk^fi weakly* in L°°(0,T;M^). 

Moreover, for any ip G Cc(]R^), there holds 

^[ fk^dv = {QifkJk),^) on (0,T), 
at j^N 

with {Q{fk, fk),f) bounded in L°°(0,T). From Ascoli's Theorem, we get 



fk'fdv 



/ (pdn{v) uniformly on [0,T]. 



As a consequence, for any given function Xe ^ C'c(IR'^ x 1^^) such that < Xe ^ 1 
and Xe{v,v^) = 1 for every (f,f*) such that |f| < e^^ and < e^^ we may pass 
to the limit (using the continuity of A = A{£,u) which is uniform on the compact 
set determined by [£^oo,^^o] and the support of Xe) 



(5.3) / D^{fk)dt — > / / \u\^A{£oo,u)xe{v,v^)diJ,dij,^dt, 

Jo k^+oD jRiVxKJV 

where we have defined for any measure (or function) A: 



D,{X) := / A(^, u) Xe{v, V.) dX{v) dX{v,). 

From the dissipation of energy p.Sj) and the estimate from below p.l9|) . there holds 

^£{t)<-D{f) with D{f):= I \ufA{£,u)ff,dvdv., 
at Jr^xr'^ 
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which in turn imphes that t D{f{t, .)) G L^{0, oo), and then 

(5.4) / D,{h)dt< / D{fk)dt= / D{f)dt 0. 
Jo Jo Jt„^ 

Gathering ()5.3|) and ()5.4p . and letting e goes to 0, we deduce that 
/ \ufA{£ao,u) dfidfi^, = on (0,T). 

The positivity ()1.18p of A{Soo,u) then imphes that fl = cSy=ijj for some measurable 
functions w : (0,T) — > and c : (0,T) — >■ Moreover, from the conservation of 
mass and momentum ()1.32p and the bound of energy ()1.33|) we deduce that c = 1 
and w = a.e. It is then classical to deduce (by the uniqueness of the limit and the 
fact that it is independent on time) that (j5.ip holds. □ 

To conclude that this weak convergence of the distribution to the Dirac mass 
as time goes to infinity implies the convergence of the kinetic energy to {i.e., the 
kinetic energy of the Dirac mass) we have to show that no kinetic energy is escaping 
at infinify as t — > T^. To this purpose we put stronger assumptions on the collision 
rate. The first additional assumption H3 roughly speaking means that the energy 
dissipation functional is strong enough to forbid it, whereas the second additional 
assumption H4 allows to use the uniform propagation of moments of order strictly 
greater than 2 to forbid it. 

Proposition 5.2 LetTc G (0, +oo] be the time of life of the solution. Then if either 
Tc < +0O, or Tc = +00 and B satisfies additional assumptions H2-H3 or H2-H4, 

we have 

(5.5) £{t) when t T^. 

Proof of Proposition 15.21 We split the proof in three steps. 

Step 1. Assume first Tc < +oo. The claim follows from the existence proof. 

Step 2. Assume now Tc = +oo and that B satisfies assumption H3: ()1.19p - p.20|) . 
We argue by contradiction: assume that S(t) -f^ 0, that is, there exists E^o > 
such that S{t) G {Soo,Sin). Reasoning as in Proposition 15.11 we get, for a fixed 
time T > and for any sequence (t„) increasing and going to infinity, that there 
exists a subsequence {tnj and a measure ft G L°°{0,T; M2) such that the function 
fkityv) := /(tnj. +t,v) satisfies ()5.2|) and 

(5.6) f Dl{f,)dt^ f Dl{-^)dt, 

Jo Jo 

where we have defined for any measure (or function) A: 

D',iX):= [ \u\'tlji\u\)xeiv,v,)dX{v)d\{v,). 
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From the dissipation of energy (jl.8j) and the estimate from below (jl.l9j) . there holds 



(5.7) ^£it) < -D%f) with D\f) := [ 



\u\^'il:{\u\) f f^dv dv^, 



which in turn implies that t D^{f{t, .)) G -^^^(0, oo), and then 
(5.8) / Dl{f,)dt< / D\fk)dt= / D\f)dt 0. 



Gathering ()5.()|1 and ()5.8|1 . and letting e goes to 0, we deduce that D^{fi) = on 
(0,T). The positivity of implies as in Proposition 15.11 that supp /2 C {0} and 
p,Sy=o. As this limit is unique and independent on time we deduce that ()5.1|) holds. 

Now, on the one hand, taking R = there holds 



(5.9) / f\v\'dv= f\v\'dv- f \v\' dv > £oo - R' > £oo/2 

J B"^ JR^ J Br 

for any t > 0. On the other hand, for T large enough, there holds thanks to ()5.1j) 

(5.10) / fdv>]- for any t>T. 

JBr/2 2 

Remarking that on -Br/2 x i?^ there holds, thanks to ()1.20|) . 

(5.11) I^N(l«l)>%^^f^l >^ii 



2 y - 4 ' 

we may put together ()5.7|1 - ()5.11|) and we get thanks to ()5.9j) and ()5.1()j) 

^£{t) < - [ [ \v - v^l'^ i){\v - v^\) f f^dvdv^ 
at JBr/Jb-^ 

< — / f dv f * I f * P dv^: < — — 

4: Jr Jrc 4 2 2 

for any t > T. This implies that S becomes negative in finite time and we get a 
contradiction. 

Step 3. Finally, assume that Tc = +oo and B satisfies assumption H4. On the one 
hand, thanks to there holds 



sup / f{t, v) \vf dv < oo. 

[0,oo) JrN 



[0,oo) 

On the other hand, arguing as in Step 2, we obtain (keeping the same notations) 
that ()5.2j) and then (from the uniform bound in Lg) 

£ifk) ^B = Sift) and Dip.) = 0. 
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The dissipation of energy vanishing imphes that 

|up/i/2* = or A(£,u) is not positive on (0, T) x M^^. 

In the first case we deduce that p, = 6^=0 as in Step 2 and then S = S(6y=o) = 0. 
In the second case we deduce, from p.lSj) . that S is not positive. In both case, 
there exists such that —>■ oo and S{Tk) —>■ and therefore ()5.2p holds since S is 
decreasing. □ 

Now we turn to some criterions for the cooling process to occur or not in finite 
time. 

Proposition 5.3 Assume that a is bounded near S = 0, and jg converges to as 
£ ^ uniformly near £ = {], then Tc = +oo. 

Proof of Proposition \5.3i It is enough to remark that, thanks to the hypothesis 
made on a and j^, the a priori bound in Orlicz norm that one deduces from (j2.13p 
as in Corollarv 12 . 71 extends to all times: 

Vt >0 WMIla < ll/inllM exp(c||/in|Lit) 

for some constant C depending on the collision rate. It shows that the energy cannot 
vanish in finite time. □ 

Proposition 5.4 Assume that B satisfies H4, that for some increasing and positive 
function Aq there holds A(£^, u) < /S.q{S) for any u E M^, £ > 0, and that /in e*" 1^''' G 
for some r > and rj G (1, 2], then Tc = +oo. 

Proof of Proposition 15.41 From the dissipation of energy ()1.8|) . the bound on A and 
the decay of the energy (|1.33p . we have 

^ > -Ao(^in) [ [ ff* dvdv, =■■ -Ao(^in) {h,R + hu) 

h,R-= I \u\^l{\u\<R} f f*dvdv^, 

h,R-= \ufl{\u\>R}ff*dvdv^. 
On the one hand, for any i? > 0, we have using ()1.32|) 



where 



h,R<R \uY f f^dvdv^ = 2R£. 

On the other hand, we infer from Proposition 13.21 (since B satisfies H4) that 

sup / f{t,v)e^''\''\" dv <Ci 
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for some r',Ci G (0, oo). Therefore 

h,R < / {A\v\^ + 4 Iv.f) 21 {\^\yR/2} f f.dvdv, 

[ {l + \vf)e'''\''\'' fdv [ {l + \v,f)f,dv,<C2e-'' 
Jrn Jrn 



< 8e 



Gathering these three estimates, we deduce 



dt 

which in turns imphes, thanks to a Gronwall argument, 

-r' R'' 

Vi?>0, inf >^i„e-^3^^-^ . 

te[o,T] R 

We conclude that £{t) > for any t G [0, T] and any fixed T > 0, choosing R large 
enough (using that ?7 > 1). O 

Proposition 5.5 Assume A(£,u) > Aq£^ with Aq > and 5 < —1/2, then Tc < 

+ 00. 

Proof of Proposition 15. 51 On the one hand, from the dissipation of energy (jl.8|) and 
the bound on A, we have 

^<-AoS' f f fU H^'dvdv,. 

On the other hand, from Jensen's inequality and the conservation of mass and 
momentum, there holds 

/ [ f f* \uf dvdv, > ( [ [ fMu\''dvdv,Y =i2£f/''. 

JrnJ^n \Jm.nJ^n J 

Gathering these two estimates, we get 

4^<-Ao^''+=^/2 
dt 

and E vanishes in finite time. □ 



Appendix: Some facts about Orlicz spaces 

The goal of this appendix is to gather some results about Orlicz spaces in order to 
make this paper as self-contained as possible. The definition and Holder's inequality 
are recalls of results which can be found in for instance. We also state and prove 
a simple formula for the differential of Orlicz norms, which is most probably not new, 
but for which we were not able to find a reference. 
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Definition 

We recall here the definition of Orlicz spaces on R-^ according to the Lebesgue 
measure. Let A : ]R_|. M.^ be a function strictly increasing, convex, such that 

(A.l) A(0) = A'(0) = 0, 

(A.2) Vt>0, A(2t) < CAA(t), 

for some constant ca > 0, and which is superlinear, in the sense that 

(A.3) 

We define the set of measurable functions / : — >■ M such that 

A{\f{v)\)dv < +00. 

Then is a Banach space for the norm 

||/|UA = inf|A>0| J^^a(W}1^ dv<l 

and it is called the Orlicz space associated with A. The proof of this last point can 
be found in [^3 Chapter III, Theorem 3]. Note that the usual Lebesgue spaces 
for 1 < p < +CXO are recovered as particular cases of this definition for A{t) = t^/p. 

Let us mention that for any / G L^(M^), a refined version of the De la Vallee- 
Poussin Theorem [23 Proposition 1. 1.1] (see also [221 1211) guarantees that there 
exists a function A satisfying all the properties above and such that 

A{\f{v)\)dv < +00. 
Holder's inequality in Orlicz spaces 

Let A be a function strictly increasing, convex satisfying the assumptions ()A.1|) . 
()A.2|) and ()A.3|) . and A* its complementary Young function, given (when A is C^) 
by 

\/y > 0, A*iy) = yiAT'iy) - A((A')-^(y)). 

It is straightforward to check that A* satisfies the same assumptions as A. Recall 
Young's inequality 

(A.4) \/x,y>0, xy < A{x) + A*{y). 

Then one can define the following norm on the Orlicz space L^*: 



Ar^*(/) = sup| / 



\fg\dv- / A{\g\)dv<l 



One can extract from Chapter III, Section 3.4, Propositions 6 and 9] the fol- 
lowing result 
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Theorem A.l (i) We have the following Holder's inequality for any f G L^, g G 



(A.5) / \fg\dv<\\fh.N^ (g). 

Jrn 

(a) There is equality in if and only if there is a constant < k* < +00 such 

that 

,A,) W.Af m V 



liA/ Viv^'Cs); \N-''-{9) 

for almost every v G M^. 
Differential of Orlicz norms 

In order to propagate bounds on Orlicz norms along the flow of the Boltzmann 
equation, we shall need a formula for the time derivative of the Orlicz norm. 

Theorem A. 2 Let A be a function strictly increasing, convex satisfying HA.l]) . 
/ TO) . / TO) , and letO < f e C\[0,T],L^) such that /(t, ■) ^ for all t G [0,T]. 
Then we have 



N^* i A' 



I/I 



-1 



(A.7) i|/,|UA 

L V VllillM. 

Proof of Theorem \A.^ From [HHl Chapter III, Proposition 6]), our assumptions on 
A imply that 

i/i 



for all 7^ / G L^. By differentiating this quantity along t we deduce: 

0. / aj,^(JIL),...^jjf,u.ljA'{^]... 



Now using the case of equality in Holder's inequality ()A.5|) we have 



9 = A' 



IM/ V VIIJIlLA. 

since the equality ()A.6|) is trivially satisfled with 

/IIm, 

and k* = N^'i^g), using that 

xy = A(x) + A*{y) 

as soon as ?/ = A'(x). This concludes the proof. □ 
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